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ABSTRACT 

We present results from spectral and temporal analyses of Suzaku and RXTE observations of the high-mass X- 
ray binary LMC X-4. Using the full 1 3 years of available RXTE/all-sky monitor data, we apply the ANOVA and 
Lomb Normalized Periodogram methods to obtain an improved superorbital period measurement of 30.32 ± 
0.04 days. The phase-averaged X-ray spectra from Suzaku observations during the high state of the superorbital 
period can be modeled in the 0.6-50 keV band as the combination of a power law with F ^ 0.6 and a high- 
energy cutoff at ~25 keV, a blackbody with ^7bb ^0.18 keV, and emission lines from Fe Kq,, O VIII, and 
Ne IX (x Lya). Assuming a distance of 50 kpc, the source has luminosity Lx ^ 3 x 10^** erg s"' in the 2-50 
keV band, and the luminosity of the soft (blackbody) component is Lbb '^ L5 x 10^^ erg s~'. The energy- 
resolved pulse profiles show single-peaked soft (0.5-1 keV) and hard (6-10 keV) pulses but a more complex 
pattern of medium (2-10 keV) pulses; cross-correlation of the hard with the soft pulses shows a phase shift 
that varies between observations. We interpret these results in terms of a picture in which a precessing disk 
reprocesses the hard X-rays and produces the observed soft spectral component, as has been suggested for the 
similar sources Her X-1 and SMC X-1. 

Subject headings: accretion, accretion disks — pulsars: individual (LMC X-4) — stars: neutron — X-rays: 
binaries 



L INTRODUCTION 

An X-ray binary system consists of a normal star orbiting 
a compact object (neutron star or black hole) to which the 
normal star transfers matter either by overflowing its Roche 
lobe or via a stellar wind. X-ray pulsars are a class of X-ray 
binaries in which the neutron star has a strong magnetic field, 
so that accreting matter follows the field lines and falls into the 
magnetic poles, generating pulses as the magnetic poles rotate 
in and out of our line of sight. While the general picture of the 
accretion mechanism is well known, the physics of accretion 
near the magnetosphere, where the neutron star's magnetic 
field begins to dominate the flow, is not fully understood. By 
studying individual X-ray pulsars in detail, we hope to gain a 
better understanding of the accretion mechanism. 

In this paper, we study the X-ray pulsar LMC X-4, which 
is a high-mass X-ray binary system consisting of a 1 .25 M© 
neutron star accreti ng from a 14.5 M© 08 III companion 
(iKelley et al.lll983t Ivan der Meer et al.l 12007 ). The neutron 
star rotates with a period of ^13.5 s, and orbit s its com- 
panion with a period of ^1.4 days (IWhitel[r978h . In addi- 
tion to the pulse period and the orbital period, the system 
has b een observed to h ave a long-term period of 30.5 ± 0.5 
days (Langet al.ni981'). Between the high and low states of 
this long-term "superorbital" period (believed to be caused 
by a precessing accretion disk that periodically obscures the 
neutron star), the X-ray flux varie s by 2 orders of mag nitude 
jHeemskerk & van ParadiisI IT9891 iNaik & Paull 120031) . The 
spectrum of LMC X-4 in the 0.1-100 keV band is weU de- 
scribed by a hard power law with a high-energy cutoff, a soft 
X-ray excess, and a strong iron emission line (INaik & Paull 
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12003^, 'La Barbera et alJ|2doTb- 

Hickox et al. (2004) showed that a soft X-ray excess was a 
common, and possibly ubiquitous, feature of X-ray pulsars, 
and that in high-luminosity systems (with Lx > lO-'^ erg s~'), 
the soft excess originates from reprocessing of hard X-rays 
from the neutron star by the inner accretion disk (Figure [l]i- 
This interpretation had been put forward previously to ex- 
plain the emission from several X-ray pulsars; for example, 
Burderi et al. (2000) suggested that in Cen X-3, a soft excess 
below 1 keV originates from reprocessing of the primary ra- 
diation by an opaque shell located at the magnetosphere. In 
this picture, the soft X-rays may be expected to show pul- 
sations that differ in shape and phase from the hard X-ray 
pulses. Endo et al. (2000) showed that the precessing accre- 
tion disk in the X-ray pulsar Her X-1 reprocessed hard X- 
rays (> 2 ke V) from the pu lsar and emitted soft X-rays (< 1 
keV), and lZane et al.l (120041) showed that the phase difference 
between hard and soft pulses changed as the accretion disk 
precesses. For the X-ray p ulsar SMC X-I, a h igh-mass X-ray 
binary similar to LMC X-4. lNeilsen et al.l (1200 4) suggested an 
explanation for the observed characteristics of the soft pulses 
in term s of precession of the accretion disk. Hicko x & Vrtilekl 
(120051) used a model of a twisted inner disk, illuminated by 
the rotating X-ray pulsar beam, to simulate pulsations in the 
soft component due to this reprocessing. They found that for 
some disk and beam geometries, precession of an illuminated 
accretion disk can roughly reproduce the soft pulse profiles. 
However, modeling the disk geometry of SMC X-1 is difficult 
since its superorbital period changes with time. The superor- 
bital period of LMC X-4 is more stable and hence may allow 
a better test of this picture. 

In this paper, we study the X-ray spectra and pulse proper- 
ties of LMC X-4 at three different phases in the high state of 
its superorbital cycle. First, we determine an improved value 
for the superorbital period based on 13 years of Rossi X-Ray 
Timing Explorer (RXTE) all-sky monitor (ASM) data and use 
it to accurately determine the superorbital phase of our three 



HUNG, HICKOX, BOROSON, & VRTILEK 



reprocessed * 
Y-^ (soft) X-rays ' 




hard pulsar beam 



illuminated disk 
region 



magnetosphere 




Fig. 1. — Schema tic of disk reprocessing in X-ray pulsars 
jHickox & VrtilelJ 120051) . Tlie hard (power-law spectrum) X-ray beam 
sweeps around and illuminates the accretion disk, which re-radiates a 
soft X-ray (blackbody) component. The observer sees pulses from both 
power-law and blackbody components, but these differ in pulse shape and 
phase. 

Suzaku observations. We then analyze the phase-averaged X- 
ray spectra and energy-resolved pulse profiles for the Suzaku 
observations and interpret them in terms of a simple model 
based on the reprocessing of hard X-rays by the precessing 
accretion disk. In Section |2] we describe the observations, 
and in Section [3] we present our analysis and results. In Sec- 
tion m we discuss the implications of our results and future 
work. All quoted uncertainties are 90% confidence for a sin- 
gle interesting parameter. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Suzaku 

Suzaku is an X -ray observatory launched in 2005 
(iMitsuda et al.ll2007n . covering the 0.2-600 keV energy range 
with two sets of instruments: X-ray CCDs (X-ray Imaging 
Spectrometer: XIS) covering 0.2-12 keV and the Hard X-ray 
Detector (HXD), a photon-counting (non-imaging) detector 
covering 10-600 keV. The HXD is comprised of two instru- 
ments, the PIN (sensitive to energies 10-70 keV) and GSO 
(sensitive to 40-600 keV). Here, our analysis focuses on XIS 
observations in the 0.5-10 keV band and the PIN in the 10-50 
keV band. We do not make use of GSO data as the signal in 
the GSO is very weak (< 1% of background). Our analysis 
uses XIS and PIN event files that have been reprocessed with 
the FTOOL aepipeline using the latest Suzaku calibration 
(version 2.2.7. 18) and screened using the standard criteria. 

Suzaku performed three observations of LMC X-4 during 
the high state of the superorbital cycle, when the source is 
bright enough to obtain a high signal-to-noise ratio. Details 
of the observations are shown in Table[T] These observations 
were carried out at the "XIS nominal" pointing position for 
effective exposures of ~20 ks. Table[T]also gives the approx- 
imate orbital phases co vered by the observ ations, estimated 
from the ephemeris of 'L evine et al.l (120001) . The exposures 
all avoid orbital phases close to eclipse (the range of (/)orb is 
0.1-0.8) and do not show any significant variations in flux 
other than the observed pulsations, so we do not expect any 
changes in absorption (from the companion star or accretion 
flow) within each exposure. 

2.1.1. XIS 

XIS consists of four detectors, each with 1024 x 1024 pixel 
X-ray-sensi tive CCDs at the foci of each of the four X-ray 
Telescopes (iNaik et al.ll2008l) . However, at the time when the 
data were taken, only three of the XIS detectors (XISO, XISl, 



Fig. 2. — Source and background regions for analysis of Suzaku observa- 
tions of LMC X-4. This image shows chip XIO for observation 1. To subtract 
background, we extracted spectra from the source region (solid) and back- 
ground regions (dashed). 

and XIS3) worked normally. Of these, XISO and XIS3 are 
front-illuminated (FI) chips, while XIS 1 is a back-illuminated 
chip. XISl has a significantly different spectral response 
from the FI chips, with greater sensitivity at energies below 
~1.5 keV and poorer sensitivity above ^ 3 keV. All XIS 
chips have an energy resolution ^ 130 eV (FWHM) at 6 keV 
(Mitsudaet al. 2007). The XIS was operated with the "1/8 
Window" option, which gave a time resolution of 1 s, cover- 
ing a field of view of 17'.8 x 2'.2. 

We used ds9 V5.6.3, an astronomical imaging and data vi- 
sualization application developed by SAO, XSELECT V2.4a, 
and XSPEC V12.6.0 for data analysis. We used XSELECT 
to combine XIS exposures taken in 3 x 3 and 5x5 event 
modes to make single event files. Next, we defined the re- 
gions for the source and the background for each observation, 
as shown in Figure |2l The source region was a rectangle of 
size 300" X 140", centered on the source and tilted 30° (for 
XISO and XIS3) and 119° (for XISl). Background events 
were extracted from two rectangles of size 260" x 130" on ei- 
ther side of the source and aligned with the source region. We 
note that the XIS background is extremely low, so that back- 
ground contributes only KiQ.5% of the total 0.5-10 keV flux 
from the source. 

For spectral analysis, we filtered the XIS events using a 
good time interval (GTI) created by merging GTIs for the 
cleaned XIS and PIN event files as w ell as the PIN back- 
ground event file (described in Section I2.1.2l i. The merged 
clean exposures range from «18 to 21 ks (Table[T]). Using the 
filtered XIS event files, we extracted spectra from the source 
and background regions and grouped the spectra so that the 
total number of events in each energy interval was roughly 
equal. The energy bands below 2.5 keV, 2.5-7 keV, and above 
7 keV are grouped by factors of 8, 16, and 32, respectively. 
Each (ungrouped) channel has a width of 3.65 eV. 

We used the tools xisrmfgen and xissimarfgen 
available in the HEASOFT version 6.9 to generate XIS redis- 
tribution matrix files and ancillary response files. Following 
the recommendations in the Suzaku ABC guid^l, we gener- 
ated 400,000 simulated photons in creating each ancillary re- 
sponse. For each observation, we combined the data from the 
two FI chips (XISO and XIS3) into a single spectrum. 

We note that "wobble" of the Suzaku spacecraft, owing to 
time-dependent thermal distortions of its chassis, can cause 
the pointing position to vary by as much as ^40". As this is a 



^ Suzaku ABC guide, 
2007-04. 



2, Section 6.6.4; JX-ISAS-SUZAKU-MEMO- 
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TABLE 1 

Suzaku Observations 








Obs ObsID 


Start Date 


'Aorb" 


XIS 


Clean Exposure 
PIN 


(ks)'' 

Merge"^ 


1 702038010 

2 702037010 

3 702036020 


2008 Jan 15 08:39:26 
2008 Feb 11 00:42:25 
2008 Apr 5 13:15:26 


0.46-0.78 
0.42-0.80 
0.11-0.64 


21.7 
23.4 

25.3 


22.9 
20.0 

24.5 


20.5 
17.9 
21.4 


■■ Based on the orbital ephemeris oflLevine et aU (2000). 











*" Note that because of gaps in the exposures, the clean exposure times are roughly half the total duration of each 
observation as indicated by the range of orbital phases. 

'^ Total good exposure time used for spectral analysis, obtained by merging the good time intervals of the XIS, PIN, and 
PIN background event files. 



significant fraction of the size of the 1/8 Window (our extrac- 
tion regions are 140" wide), wobble may cause source flux 
to be lost outside the window during the XIS exposures. We 
examine the thermal wobble by deriving the centroid position 
of the source in detector coordinates (DETX and DETY in the 
event files) in bins of size 50 s. For > 97% of the total ex- 
posure time the thermal wobble is <20", for which there is 
negligible flux loss. Observation 3 has a brief period (<2 ks) 
of wobble >30". We derive a rough estimate of the flux lost 
outside the window during these periods, by calculating the 
observed source counts after shifting the extraction region by 
an amount and direction corresponding to the observed ther- 
mal wobble. At most only 15% of the flux is lost outside the 
window at any given time, and the total flux lost over the full 
observation is <1.5%. We thus verify that losses due to ther- 
mal wobble are minimal, and for simplicity we do not exclude 
the brief intervals of large wobble from the analysis. 

2.1.2. PIN 

We include PIN data in the spectral and timing analysis. 
As the HXD is a non-imaging instrument, we cannot esti- 
mate background directly from the observations. The PIN 
background consists of an instrumental component as well as 
the cosmic X-ray background (CXB). The instrumental back- 
ground is estimated using the version 2 "tuned" background 
event filefl Spectra were extracted for both the observations 
and the background file and corrected for dead time. As with 
the XIS, spectra were extracted using the merged GTI de- 
scribed above. The CXB component is modeled during spec- 
tral analysis, as a power law with a high-energy exponential 
cutofQ this component is only ~2% of the observed flux. 

2.2. RXTE 



RXTE was launched on 1995 December 30 ("Swank 2006^. 
The RXTE ASM observes bright X-ray sources in the range 
of 1.5 - 12 keV, to explore the variability of sources including 
black holes, neutron stars. X-ray pulsars, and bursts of X-rays. 
ASM scans about 80% of the sky every orbit, allowing mon- 
itoring on timescales of 90 min utes or longer . Details about 
the RXTE ASM can be found in lLevineet aH ([199 6). In this 
paper, we use 13 years of the RXTE ASM data on LMC X- 
4 from 1996 to 2009. We retrieve the one-day average light 
curve from the public MIT Web sitJH. 

3. DATA ANALYSIS AND RESULTS 



« JX-ISAS-SUZAKU-MEMO-2007-04 

^ ftp://legacy.gsfc.nasa.gOv/suzaku/data/background/pinnxb_ver2.0/ 

'" Suzaku ABC guide, v. 2. Section 7.3.3.1 

'^ http://xte.mit.edu/ASM_lc.html 



3.1. Superorbital Period Determination 

In order to determine the superorbital period during the time 
of our Suzaku observations, we first apply an epoch folding 
method to the 13 years of RXTE data to obtain an updated 
superorbital period. By comparing the variance of the folded 
light curve with the varian ce within each bin (the ANOVA 
method; D avieslll990lll99ll) . we calculate the Davies L statis- 
tic for different superorbital periods. Gaussian fits to the pe- 
riod versus L statistic, for different numbers of phase bins, 
all give about 30.316 days at high significance. A power 
spectrum using the Lomb Normalized Periodogram (IScargld 
119821) is also consistent with 30.32 days. This result is con- 
sistent with the supero rbital period determined using five 
years of RXTE data by iPaul & Kitamotol ( 120021) . This sug- 
gests no significant change in the long-term period over the 
13 years of RXTE obse rvations, in contrast to the results of 
iPaul & Kitamotol (l2002h who had found a long-term decrease 
in the period over the preceding 20 years. We estimate the er- 
ror by re-folding the light curve into various numbers of bins 
and examining the resulting best fit. Both the ANOVA and 
Lomb Normalized Periodogram methods determine a long- 
term period of 30.32 ± 0.04 days. Figure |3] shows a plot of 
our best-fit period superposed on the 1 3 year RXTE ASM light 
curve. A study was made of possible variability in the long- 
term period over the RXTE ASM lifetime and the folded light 
curve gave similar results near the times of the Suzaku obser- 
vations as the steady period model. 

We note that the L statistic does not make use of uncertain- 
ties in the RXTE lightcurve. As a check, we performed an 
e-folding analysis using the FTOOL ef search to derive the 
period that maximizes x^ relative to no variation. We used 
8, 16, or 32 phase bins with time resolutions of 0.09, 0.05, 
or 0.02 days, respectively. We obtained long-term periods of 
30.29 ±0.09, 30.29 ±0.05, and 30.32 ±0.02 (where the er- 
rors are the period resolutions of the search), corresponding to 
maximum x^ of 890, 959, and 1070, respectively. These esti- 
mates are consistent with our result derived from the ANOVA 
and Lomb Normalized Periodogram methods. 

Our improved fit to the superorbital period allows us to ac- 
curately determine the phases of our Suzaku observations as 
shown in Figure |4] We define the start of the high state of the 
superorbital cycle (03o = 0) to be at JD 2454560.0. 

3.2. Phase-averaged Spectroscopy 

For each observation, we perform spectroscopic fits to the 
XIS FI, XIS 1, and PIN spectra simultaneously, for a com- 
bined spectral range of 0.6-50 keV. We fit all three spectra 
with the same spectral model, but allow the overall normaliza- 
tions to vary in order to allow for uncertainties or systematic 
offsets in the flux caHbration. Fits to the power-law spectrum 
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Fig. 3. — One-day average light curve of LMC X-4 covering 13 years of observations (1996-2009) by the RXTE ASM. The magenta line is our best fit to the 
30.32 it 0.04 days superorbital period of LMC X-4. The hats on the error bars denote the widths of the time over which the data were averaged (1 day for the 
ASM data). 



of the Crab nebula indicate that the best-fit normalizations for 
the three XIS chips are consistent to within a few percent, 
while the best-fit normalization for the PIN relative to XIS is 

i.i(Q. 

We first fit the phase-averaged spectra with a simple contin- 
uum model of power-law plus blackbody that ha s previousl y 
been successfully employed for LMC X-4 (Pau let alJl2002l) 
as we ll as a number of other X-ray pulsars (IHickox et alJ 
12004ft . Previous studies have found evid ence for a high- 
energy cutoff in the power-law component; iLa Barbera et alJ 
I use a cutoff of the form 



l/(l+e' 



,(£'-£'(;ul)/£'fO 



'), 



(2) 



„(E-E, 



'cui)/£roid ■£'>£' 



(1) 



with Ecut ~ EfoU ~ 1 8 keV. We find that for the Suzaku data 
the spectra are better described by a F ermi-Dirac cutoff (e.g., 
ICoburn et alJl2(J()2LlRivers et alJl20Toh of the form: 



JX-ISAS-SUZAKU-MEMO-2008-06 



with Ecut ~ 27 keV and Zifoid ~ 8 keV. 

In the XIS data, we exclude the energy range 1.5-3 keV 
to avoid in strumental effects near the Au and Si edges 
jMiller et al . 2008; Sae z et a l. 2009). We include neutral hy- 
drogen absorption fixed to the Galactic value of A^h = 0.057 x 
10^^ cm"^ (.Paul et al] 120021) . Our best-fit parameters for 
this simple model give a blackbody with kT^B ^0.18 keV 
and a power law with pho ton index F ^ 0. 6, similar to re- 
sults f rom previous studies dPaul et a l. 2002; La Barbera et alj 
120011) . Spectral fits are shown in Figure |5] and fit parameters 
are given in Table |2] 

A single cutoff power law provides a reasonable fit to the 
high-energy component in the full range from 1 to 50 keV; 
the best-fit photon index for the joint XIS and PIN analysis is 
consistent (within 2%) from that derived from fitting the XIS 
data alone. We note that Figure |5] shows a departure from the 
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Fig. 4. — Close up of ASM count and the region of the Suzaku observations. Blue arrows indicate the times of our Suzaku observations. The superorbital 
phases for the tliree observations (left to right) are 03o = 0.39, ^30 = 0.27, and 4>3Q = 0.07. The zero superorbital phase is defined as the beginning of a high state 
(JD 2454560.0). The magenta line is our best fit of the superorbital period of LMC X-4. The hats on the error bars denote the widths of the time over which the 
data were averaged (1 day for the ASM data). The widths of the arrows denote the length of our observations (about 20ks each). 



cutoff power law at ^^40 keV for observation 1 . This feature 
does not appear in observations 2 and 3, and so for simplicity 
we model all three observations with the simple cutoff power 
law. 

In the XIS spectrum at £ < 10 keV, it is clear from the 
residuals in Figure |5] that other, line-like spectral features are 
present in addition to the simple continuum model. Emission 
line features are often observed in X-ray binaries; for exam- 
ple, in SMC X-1, several emission lines from O, Ne, Mg, Si, 
S, Ar , and Fe were needed to fit the spectrum ( Vrtilek et alj 
12005b . In addition, using the Chandra High-Energy Trans- 
mission Grating Spectromet er and the XMM-New ton Reflec- 
tion Grating Spectrometer, iNeilsen et alJ (l2009h performed 
X-ray spectroscopy of LMC X-4 and detected several X-ray 
emission lines. In the Suzaku observations, we find at least 
four emission lines including Ovill, NelX (x Ly^), Fe K^, 
and broad Fe (line energies and equivalent widths (EWs) are 
given in Table O. The total x^ values for every observation 
are lower by at leas t 30 when a ddin g each individ ual emis- 
sion line. Following lAvnil (1 1 97 6") and'Y aqoobI (I1998D . we are 
>90% confident that the line features are detected. The av- 
erage x^ per degree of freedom for the three observations is 
reduced from 2.5 to 1.2. The spectral fits are shown in Fig- 
ure |6] and the results are listed in Table ID 

We note that for observations 1 and 3, with the inclusion 
of emission lines, the best-fit relative normalizations between 
the PIN and the XIS FI spectra are broadly consistent with the 
value of 1 . 16 derived from calibration studies and higher than 
the value of « 1 .07 for the pure continuum model. 

As an alternative to blackbody model for the soft com- 
ponent, we also consider models for optically thin thermal 
emission, which have previously been used to describe the 
soft e mission from LMC X-4 (Paul et al. 2002; Naik & Paul 
12004ft . We consider both a thermal bremsstrahlung model and 
an optically thin (MEK AL) model that include s contributions 
from metal lines (e.g.. iHaberl & PietschI 120051) . We perform 
fits using a model identical to that described above including 
emission lines, but replacing the blackbody component with 
the thermal bremsstrahlung and MEKAL models. We find 
that these can represent the soft component equally well to 
the blackbody model (although the MEKAL model requires 



very low metallicity < 2% of solar). Therefore, we cannot 
rule out these models for the soft excess on spectral grounds. 
However, we can rule them out on physical grounds, follow- 
ing the argument of (Hicko x et a l. 2004). It is p ossible that the 
diffus e cloud of gas around the neutron star dBoroson et alj 
1200 Ih would reprocess the hard X-ray and produce optically 
thin emission that might be described by optically thin emis- 
sion models. However, in order to produce high Lsott as ob- 
served in LMC X-4, the system must have the emission region 
^cioud > 10'^ cm, which is larger than the size of the emission 
region /?cioud ^10" cm as inferred from the decrease in soft 
flux during the eclipse. Further, such a large region cannot 
produce soft pulsations as observed. Therefore, although the 
optically thin models provide an adequate description of the 
data, we prefer the blackbody model for our physical interpre- 
tation of the soft component. 

For all the spectral fits shown in Figs. |5] and |6] the black- 
body component dominates the energy range below 1 keV, 
while power-law emission dominates above 2 keV. We pro- 
pose that the hard (power-law) X-rays from the pulsar are re- 
processed by the inner accretion disk, where the soft pulses 
are emitted. The power-law flux is (8-10) x 10"^ erg cm~^ 
s"', assuming the absolute flux calibration of the XIS FI 
chips. Taking the distance to LMC to be 50 kpc, we obtain 
Lx ^ 3 X lO-'^ erg s"' from 2 to 50 keV. The bolometric lumi- 
nosity of the blackbody component (taken directly from the 
XSPEC normaHzation) is ^ 1.5 x 10^^ erg s"'. 

Using the observed temperatures and luminosities, we can 
derive the size of the blackbody-emitting region, assuming 
that the emission is from a spherical region that partially cov- 
ers the central source. Note that this differs from another com- 
monly used estimate for the blackbody radius, which assumes 
a circular emission geometry. We prefer the spherical par- 
tial covering geometry, because for this source the blackbody 
emission is likely due to the reprocessing of the hard emission 
from the neutron star (Hickox et al. 2 004). For a partial cover- 
ing spherical shell and an incident hard X-ray l uminosity Ly, 
the ra dius of the shell (as shown in Figure 9 of iHickox et alJ 
llOOl is 



^BB = 



AnaT^^- 



(3) 
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TABLE 2 
Fits Using a Power Law and a Blackbody to Suzaku Observations 



Parameter 


Observation 1 


Observation 2 


Observation 3 


Xl (317 dof) 

Wh(x 10^2 cm-2) (fixed) 


2.90 

5.7 X 10-- 


2.53 
5.7 X 10-2 


2.18 
5.7 X 10-2 


Continuum 



kTsB (keV) 

A'bb" 
Power-law F 

£cut (keV) 
fifoid (keV) 



0.181 ±0.002 

+0.09 1 
-0.08-' 

0.56 ±0.02 

(0.67 ± 0.02) X 10-2 

25.3+'o 

,+0.6 



7.1: 



-0.5 



0.174 ±0.002 
(6.99±0.11)x lO-"* 

0.56 ±0.02 

(1.05±0.03)xl0-2 

24.8^1;3 

8.6 ±0.6 



0.181 ±0.002 

(5.57+«»?) X 10-4 

0.56 ±0.02 

(0.83 ± 0.02) X 10-2 

24.7!|-2 

8.1 ±0.6 



Relative Normalizations 



Axisi 

>lpiN 



1.01 ±0.01 
1.07 ±0.03 



0.95 ±0.01 
1.05 ±0.03 



1.02 ±0.01 
1.09 ±0.03 



' XSPEC normalization, units of 8.36 X 10 * erg cm 2 



s (bolometric). 



' XSPEC normalization, units of photons cm s at 1 keV. 



TABLE 3 
X-RAY Emission Lines Used in Fitting Spectra of LMC X-4 



Line 



Wavelength (A) 



Line 
Energy (keV) 



Observed 
Energy (keV)" 



Equivalent 
Width (eV) '' 



O VIII Lye 
NelX (NexLyc) 
FeKa 
Broad Fe 



18.97 

13.45(12.13) 

1.94 

1.89 



0.65 
31 (1.02) 
6.38 



0.60-0.62 
0.98-0.99 
6.40-6.51 
6.52-6.61 



41 ±13 

13±5 

25 ±19 

150 ±50 



" Range of best-fit line energies for the three observations. 

'' Average EWs and uncertainties between the three observations. 

Using Lx ~ 3 X lO"*** erg s"' and kT^B '^ 0.18 keV, we have 
Rbb ~ 1.5 X 10^ cm. In this picture, the covering fraction O 
of the spherical shell is given by simply fl = Lbb/Lx and so is 
relatively small (fl ~ 0.05). 

In LMC X-4, an estimate for the magnetic field 
comes from the observed cyclotron feature at ^ 100 
keV (iLa Barbera et all 1200 Ih . We note that the exis- 
tence of this feature has been a matter of debate (e.g., 
IXsygankov & Lutovinov 2005), but since it provides the only 
existing estimate for the magnetic field in this source, we 
adopt it in the subsequent analysis. A cyclotron energy of 100 
keV would indicate B ^ lO'^ G. An estimate of the magneto- 
spheric radius (e.g.. lGhosh & Lambl[T97l iFrank et al.ll2002l) 
is 



Ru 



' 0.5Ra 



1.5xlO'^m;/X°^VXrcm, (4) 



where Ra is the standard Alfven radius. Here, nii is the mass 
of the neutron star in Mq, Rf, is its radius in 10'' cm, L37 is 
the X-ray luminosity in 10"^^ erg s"', and B12 is the neutron 
star surface magnetic field in 10'^ G. Assuming mi = 1.25, 
R(, r^ 1, L37 = 30, and Bn ~ 10, we obtain R„, '- 2 x 10^ cm 
and Rbb /R,,,-^ 0.1. 

We note that if we evaluate Lx only in the range 0.5-10 
keV, we obtain Lx '^ I x 10^^ erg s"' and Rbb ^^9 x 10^ cm, 
similar to r esults obtained by ASCA in the same energy range 
JPaul et al.| f2002): this highlights the need to consider higher 
energies where the bulk of the hard component is emitted. 

3.3. Tuning Analysis 

Our pulse-phase-averaged spectral analysis suggests that 
the soft and hard spectral components for LMC X-4 have dif- 



ferent physical origins (as discussed in detail bv lHickox et al.l 
120041) . We examine the time variation of these components 
separately, by deriving energy-resolved pulse profiles at the 
different superorbital phases. We derive pulse periods from 
the PIN observations, which have time resolution of 61 /is. 
We are also able to study pulsations with XIS, which in 1/8 
Window mode has time resolution (Is) that is sufficient to 
produce pulse profiles for LMC X-4's 13.5 s period (the XIS 
event arrival times were randomized within the 1 s time bins). 
Any variations in source flux on timescales longer than the 
pulse period should not significantly affect the derived pulse 
profiles. Therefore, to maximize signal to noise in the timing 
analysis, we used the full clean GTIs for each detector (unlike 
for the spectral analysis, for which we filtered the XIS and 
PIN data with a common GTI). 

We corrected both the PIN and XIS events to the rest 
frame of the solar system using the aebarycen tool and 
corrected for the orb ital motion o f the LMC X-4 system us- 
ing the ephemeris of iLevine et al.l (^00). Using the FTOOL 
efsearch, we derived pulse periods for the PIN data by 
searching for the maximum in the x^ versus folding period. 
We obtain Ppuise = 13.5087 ±0.0002 s, 13.5091 ±0.0001 s, 
and 13.5096 ±0.0001 s for observations 1-3, respectively. 
Uncertainties are estimated both from formal errors on the 
fit to the x^ versus period and also by a jackknife resampling 
technique, calculating the period after excluding successive 
^4 ks GTIs; these two estimates give consistent uncertain- 
ties. The drift in Ppuise over the three months between obser- 
vations 1 and 3 is consistent with typical ra tes of spin evolu- 
tion observed prev iously for LMC X-4 (e.g.. lWoo et al.llT996t 
ILevine et al.ll2000l) . 
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TABLE 4 
Fits Including Emission Lines to Suzaku Observations 



Paiameter 


Observation 1 


Observation 2 


Observation 3 


Xi (308 doO 

A'h (X 10^2 cm-2) (fixed) 


L27 
5.7 X 10-2 


L20 
5.7 X 10-2 


1.08 
5.7 X 10-2 


Continuum 



^rBB (keV) 
Power-law F 

£cut (keV) 
£fold (keV) 



+0.003 



"•1^3-0.002 

(5.69!°|2) X 10-^ 

0.64 ±0.02 

(0.73 ± 0.02) X 10-2 

26 7+*' 
''"■'-i.o 



6. 



)+0.6 
'-0.5 



0.174 ±0.003 

+0.16-, 
-0.1 3-' 

0.63 ±0.02 

(1.12±0.03)x 10-2 

26.5t\i 

8.3 ±0.6 



0.183 ±0.003 

(5.05 ±0. 12) X 10-* 

0.62 ±0.02 

(0.88:003) X 10-2 

26.f!j:J 

7.9 ±0.6 



Emission Lines 



0.62 ±0.02 
5.0 X 10-2 

(2.5!«-^) X 10-3 
0.98 ±0.02 
5.0 X 10-2 

(3.8!g-^) X 10-* 
6.41 ±0.03 
5.0 X 10-2 

(7.0!2f ) X lO-'^ 

6 53+00'' 

-fiQ5 

37+008 

(3.7!^j) X lO-'' 



0.60 ±0.02 

5.0 X 10-2 

(3.4;i|) X 10-3 

99+0.02 

5.0 X 10-2 

(4.7:0 ») X 10-* 

6.40 ±0.04 

5.0 X 10-2 

(7.9:^'5) ^ 10-5 
6 61+^07 

37+''-0' 
(4.5:0;9) X 10^ 



0.61 ±0.02 

5.0 X 10-2 

(2.5!<!-^) X 10-3 

0.98 ±0.03 

5.0 X 10-2 

(2.4 ± 0.7) X 10-" 

6 SI +006 

5.0 X 10-2 

(6.5!2-3) X 10-5 

6.52 ±0.09 

53+*- '5 
(3.7!<!J) X 10-" 



£o (keV) 

fTo (keV) 

/q (photons cm-2 s-' ) 

£Ne (keV) 

o-Nc (keV) 

^Ne (photons cm-2 ^-1 ^ 

iipe (keV) first Gaussian 

(TFe (keV) 

/pc (photons cm-2 ^-1^ 

iiPe (keV) second Gaussian 

crp, (keV) 

/pc (photons cm-2 ^-1^ 



Relative Normalizations 



Axisi 

ApiN 



1.01 ±0.01 



1.15 



+0.04 



0.95 ±0.01 
1.12 ±0.03 



1.01 ±0.01 
1.15 ±0.03 



" XSPEC normalization, units of 8.36 X 10 ° erg cm 2 s ' (bolometric). 
'' XSPEC normalization, units of photons cm-2 g-l ^j j j^gy 

As a check, we derived pulse periods from the XIS data, 
using the combined event file for the XISO, XISl, and XIS3. 
We obtain similar pulse periods of Ppuise = 13.5089 ± 0.0001 
s, 13.5095 ±0.0001 s, and 13.5100 ±0.0003 s; the resulting 
pulse profiles are qualitatively identical if we use these peri- 
ods rather than those from the PIN. 

We note that we have not ac counted for any d eviation 
from orbital ephemeris derived bv lLevine et alj (120001) . which 
could possibly affect the derived spin periods. However, the 
relative pulse profiles in the different X-ray bands, which are 
the focus of this study, will be largely independent of the ex- 
act absolute spin period. We leave a detailed analysis of the 
spin evolution of LMC X-4 to a future study. 

Using the above spin periods, we created separate pulse 
profiles for the XIS in the 0.5-1 keV ("soft") band (which 
is dominated by the soft blackbody spectral component) and 
the 2-10 keV ("medium") band (which is dominated by the 
hard power-law component). We also derived profiles for the 
PIN in the 10-50 keV ("hard") band that is entirely domi- 
nated by the hard (cutoff) power law. We have subtracted the 
background flux in producing the PIN pulse profiles (the XIS 
background is negligible). 

The pulse profiles are shown in Figures |7]-(9] and show sig- 
nificant differences in the pulse profile for the different X-ray 
energies. The soft and hard pulses in all three observations 
are broad and show single peaks, while observations 2 and 
3 show double-peaked profile s in the medium band, similar 
to those observed in Her X-1 (Zane et al.' 2004) and SMC X- 
1 (Neilsen et al. 2004: Hickox & Vrtilek 2005) . although the 
pulsations for observation 3 are relatively weak. In contrast 



observation 1 shows only a single peak in the medium pulse 
profile, similar to the observed soft pulses. The variations be- 
tween the hard and medium pulses indicate that the spectral 
shape of the medium component may change with pulse phase 
and further that these variations may differ between observa- 
tions. 

In addition, the phase shifts between pulses at different en- 
ergies vary among the observations. In observation 1, all three 
pulses are almost in phase with each other, while observations 
2 and 3 show notable phase shift. We study this phase shift 
in more detail by cross-correlating the soft profiles against the 
hard and medium profiles (see Figure [TOb. We first examine 
the cross-correlation between soft (0.5-1 keV) and hard (10- 
50 keV) pulses, shown in the top panel of Figure [TOl These 
energy bands represent the bulk of the blackbody and hard 
power flux, respectively. Observation 1 shows a strong pos- 
itive cross-correlation between soft and hard pulses peaking 
at ^^-40°, indicating that the soft pulses slightly lag the hard 
pulses. The cross-correlations for observations 2 and 3 peak 
almost exactly out phase at ^180°, although there is an indi- 
cation that the observation 3 may be slightly shifted to shorter 
phase lag (by ^10°). However, the precise relative lag is diffi- 
cult to measure because the pulse profiles have relatively few 
(10) bins. 

The cross-correlations between the medium (2-10 keV) 
pulses and the hard and soft pulses are shown in the bot- 
tom two panels of Figure [TOl For observation 1, the cross- 
correlations peak near zero as the phase of the medium pulses 
closely aligned with the soft and hard pulses. The cross- 
correlations for observations 2 and 3 are more complex, be- 
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Fig. 5. — Phase-averaged Siizaku spectra for the three observations, fitted 
with a continuum consisting of a blackbody (dashed) plus power law with 
high-energy cutoff (dash dot). The FI, XISl, and PIN spectra are shown in 
blue, red, and green, respectively. The individual spectral components are 
shown, including the CXB model for the PIN spectrum. Note the significant 
line-like residuals. 

cause of the multiple peaks in the medium pulse profiles. We 
note that the phase offset between medium and hard pulses 
varies more strongly than between medium and soft pulses. 

4. DISCUSSION 

We use 13 years of available RXTE/ ASM data to determine 
a superorbital period of 30.32±0.04 days for LMC X-4. The 
long time series allows an improvement by a factor of 10 in 



Fig. 6. — Phase-averaged spectra from Suzakii observations with blackbody 
(dashed), power law with high-energy cutoff (dash dot), and emission line 
(dotted) components shown. The FI, XISl, and PIN spectra are shown in 
blue, red, and green, respectively. The emission lines are broad Fe, Fe Kq, 
O Vlll, and Ne IX (X Ly«). The inclusion of the lines significantly improves 
the residuals compared to the continuum fits shown in Figure|5] 



precision compared to the measurement of lLang et al.' ('1981|). 
Our period estimate is consistent with that of Lang et al] 
([1981), indicating that the period has been relatively stable 
over the past ^-^30 years. The improved period allows us to 
accurately determine the superorbital phases of our Suzaku 
observations and will facilitate future studies that require ac- 
curate estimates of (/>3o. 



SUZAKU X-RAY SPECTRA AND PULSE PROFILES OF LMC X-4 




1.1 



K 1.0 



o 



0.9 



0.8 



ri^ 



r 



^ 



r 



^ 



2-10keV(XIS) 
ilO-50 keV (PIN) 




o 



1.0 
Pulse phase 

Fig. 7. — Energy-resolved pulse profiles for observation 1. The soft (0.5- 
1 keV, red) pulses are dominated by the blackbody component, while the 
medium (2-iO keV, blue) and hard (10-50 keV, black) pulses are dominated 
by the power law. The profiles for observations 2 and 3 are shown in Fig- 
ures[8]and|9] respectively; note the variation in the medium pulse profile and 
the varying phase offset between the hard, medium, and soft pulses for the 
different superorbital phases. 

In the Suzaku observations, we find that the continuum over 
roughly two decades in energy (0.6-50 keV) can be described 
by a blackbody plus a cutoff high-energy power law. In ad- 
dition, we detect four emission features as listed in Table [3] 
(All o f the four ernission lines have been detected in SMC 
X-1; dVrtilek et all l2005h .) Limited by the energy resolu- 
tion (~ 130 eV) of XIS detectors, we cannot clearly iden- 
tify the line with energy 0.98 keV. Our possible candidates 
are Neix (0.91 keV) and Nex Ly„ (L02 keV). For the Fe 
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Fig. 8. — Energy-resolved pulse profiles (as in Figure|7), for observation 2. 

Kq, li ne, it has been detected from many previ ous studies 
(e.g., lLaBarberaetal."2001"; 'Naik & Pauj |2003l and refer- 
ences therein). Neilsen et al. (2009) suggested that Fe Kq, line 
originates at the outer accretion disk. In addition to Fe Kq, at 
6.4 keV, we also detect a broad Fe line. Using M =1 .25Mq for 
the neutron star and Rbb ^ 1 .5 x 10^ cm, we find the velocity 
at the inner accretion disk to be ~ 1.1 x 10"^ km s"', which 
corresponds to Doppler shift of 0.2 keV for Fe K^. This is 
of similar order (although somewhat smaller than) the aver- 
age <T w 0.4 keV derived from the spectral fits. Given this 
rough correspondence (and the considerable uncertainties in 
the geometry of the accretion flow), we argue that this broad 
Fe line is Doppler broadened Fe Kq emitted from the inner 
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Fig. 9. — Energy-resolved pulse profiles (as in Figure|2), for observation 3. 



accretion disk due t o the high velocity, as also suggested by 
lNeilsenetal"](l2009h . 

It is challenging to estimate the fluorescent Fe line flux that 
would be expected for reprocessing by the inner disk. For 
the partial spherical geometry described ab ove with fj = 0.05, 
a mo del of reflection from a neutral disk (iGeorge & FabianI 
119911) predicts an Fe K line EW - 25 eV, similar to the EW of 
the observed narrow line but smaller than the broad Fe feature, 
which has EW = 150 ± 50 eV. We note, however, that such a 
broad feature could be produc ed by emission from iron in a 
range of ionization states (e.g.. lRamsav et al.l2002l) and could 
also be due to complex absorption around the Fe K edge, as 
observed in some active galactic nuclei (e.g., .Yagoob et all 
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Fig. 10. — Cross-correlation between energy-resolved pulse profiles, show- 
ing the variation in relative phase and shape of the pulses for the three Suzaku 
observations. A positive phase shift coiTesponds to the first pulses listed lead- 
ing the second pulses listed. For example, in the top panel, the peak at ~ ^0° 
shift for observation 1 indicates that the soft pulses lag the hard pulses by 
~40°. 



Il995h and polars (e.g.. iDone & Magdziarzlll998l) . Finafly, 
disk reflection would predict the presence of a Compton re- 
flection feature which we have not included in our simple 
power-law model for the continuum, and which might affect 
the observed EW. In light of these uncertainties, it is notable 
that a broad Fe feature with similar EW ^ 150 eV has also 
been observed in Her X-1, for which correlations in the pulse 
profiles for the soft (blackbody) component an d the Fe K line 
suggest that both originate from the inner disk (Ramsay et alj 
I2OOZ ; , Zane et al. : 2004) . We conclude that inner disk repro- 
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Fig. 1 1 . — Model of the inner regions of a twisted, warped accretion disk, 
illuminated by the neutron star and observed at three superorbital phases cor- 
responding to the Suzaku observations. We set the inclination of the disk to be 
70° >van der Meer et al.i2007iy The orange shaded regions show areas on the 
disk which are illuminated by the central star and are visible to the observer. 



cessing is a plausible explanation for the observed Fe line 
emission, although further work is needed to robustly deter- 
mine the physical origin of the broad feature. 

The pulse profiles of LMC X-4 show that the profiles at dif- 
ferent energies have different pulse behaviors. In all three ob- 
servations, the soft (0.5-1 keV) and hard (10-50 keV) profiles 
show single-peaked, smooth, roughly sinusoidal pulsations. 
Single-peaked soft pulses have been observed in the well- 
studied X-ray binary pulsar systems, Her X-1 and SMC X-1 
( IZane et alj|2004 ' 'Hickox & Vrtileai2005l) . iHickox & Vrtileg 
(120051) using a simple pencil beam model with a twisted in- 



ner accretion disk were able to reproduce the observed single- 
peaked sinusoidal pulsations. We suggest that the soft pulses 
of LMC X-4 are also due to reprocessing of hard X-rays in 
the inner accretion disk. 

The medium (2-10 keV) component, however, shows 
single-peaked pulses in observation 1 and double-peaked 
pulses in observation 2 and (although weakly) in observa- 
tion 3. Multiple-peaked pulses at E > 2 keV have bee n re- 
ported several times in the literature; IWoo et al] (119961) ob- 
served double-peaked pulses with one dominant peak, while 
|Naik&Paul (2004) and Paul et al. (2002) find pulses with 
three equally strong peaks. Weak o r abse nt pulsations have 
also been observed (L a Barbera et alj2001 l). 

Further, the phase shift between the pulsations changes sig- 
nificantly between the three observations. In particular, in ob- 
servation 1 the pulses in all three bands are largely in phase, 
while for observations 2 and 3, the hard and soft pulses are al- 
most 180° out of phase, while the medium pulses are shifted 
somewhat less relative to the soft pulses. 

Fol l owing previous st udies of Her X-1 (e.g.,'Ramsav et 
200l IZane et alj l2004h and SMC X-1 (Hickox & Vrtile" 



20051) . we attempt to interpret the observed pulse profiles 



in terms of the soft component being reprocessed emission 
from the inner disk. We begin by considering a fiducial 
model for the inner regions of the accretion disk at the su- 
perorbital phases of the three observations (Figure [TTTi. The 
model of a warp ed, twisted disk is the same as that used 
by Hic kox & Vrt ilek (2005) to describe the superorbital vari- 
ation in SMC X-lH We assume a disk inclination of 
70°, consistent with the orbital inclination of LMC X-4 
dvan der Meer et al] 120071) . The area of the disk illuminated 
by the neutron star is shown in the shaded region. While 
we do not have strong constraints on the shape of the inner 
disk for LMC X-4, this model gives us a qualitative picture of 
what we might expect for reprocessing by a warped, precess- 
ing disk. 

Another con sideration is the g e ometr y of the pulsar beam. 
For SMC X-1. IHickox & Vrtilekj (l2005h found that pencil- or 
fan-shaped beams from the neutron star can provide an ac- 
ceptable description of the observed pulse profiles, although 
the pencil beam is preferred at some superorbital phases. 
For the super- Eddington luminosities observed for SMC X-1 
and LMC X-4, theoretical models generally predict that ra- 
diation cannot escape from the top of the accretion column 
and so favo r a fan beam emitted from the wall s of the col- 
umn (e.g. Beckei 119981: iBecker & WoIfS 120051) . However, 
fan beam models universally predict a double- or quadruple- 
peaked hard pulse profile, in contrast to the single peaks ob- 
served in LMC X-4. Further, fan beams can occupy a wide 
and complex range o f geometries (for example the reverse fan 
beam suggested by Scott. Leahy. & Wilsonll20()0l for Her X- 
1), and a full treatment requires a detailed description of the 
accretion flo w as well as effects such as gravitational light 
bending (e.g.. lMeszaros & Rifferti[r988: Leahy .2003). Given 
these difficulties, in the following discussion we consider only 
a simple pencil beam and leave a precise determination of the 
beam geometry for future studies. 

Assuming that the radiation pattern consists of two pencil- 
shaped beams pointing in opposite directions, then if the 
beams are not exactly in the plane of the accretion disk, one 
beam will produce a strong observed hard pulse (as we ob- 

'^ The disk corresponds to the best- fit model for a pencil beam, with pa- 
rameters given in Table 2 of lHickox & Vrtile k 1,2005). 
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serve for LMC X-4) while the other beam may more strongly 
illuminate the disk. In observation 1, the largest visible area 
of the disk is almost exactly in the line of sight, and we may 
expect the hard and soft pulses to be almost in phase with each 
other. Figure nn shows that as the disk precesses, the largest 
visible area of the disk is shifted to the side, which would 
suggest that the phase shift between the pulsations from the 
neutron star and those from the reprocessing disk should vary 
with superorbital phase. This is indeed what we observe for 
the shift between the medium and soft X-ray pulses (shown 
in the bottom panel of Figure [TOl l. for which the peak of the 
cross-correlation increases monotonically between the three 
observations. Such variations in the ^2-10 keV and ~0.5-l 
keV pulse profile s have been interp reted similarly for Her X-1 
dZane et alj|2004 and SMC X-1 ( Neilsenet al.]|2004b . 

However, these studies were limited to X-ray data below 
10 keV. With the existence of the higher-energy PIN observa- 
tions, it makes sense to consider the hard (10-50 keV) emis- 
sion to be the best indicator of the luminosity in the neutron 
star beam, as the spectrum peaks at these energies. As dis- 
cussed above, the hard and medium pulses show significantly 
different behavior. The cross-correlations in Figure [10] indi- 
cate that the hard pulses are either close to in phase (obser- 
vation 1) or close to completely out of phase (observations 2 
and 3). 

This behavior suggests a different interpretation for the 
variation in the phase offsets with 03o. One idea is that if the 
beam aligned with the observer also illuminates the largest 
visible area of the accretion disk, then we might expect the 
hard and soft pulses to be almost exactly out of phase. It is 
therefore possible that between observations 1 and 2, the hard 
X-ray beam and the visible region of the disk moved abruptly 
relative to each other, thus changing which of the two pulsar 
beams is illuminating the visible region of the disk. In terms 
of this picture, the observations would be consistent with a 
simple "donut" geometry for the reprocessing region, rather 
than the more complex twisted disk as shown in Figure [TT] 

We further note that the amplitudes of the soft pulses are 
equal to or slightly weaker than those of the hard pulses. If 
the pulses arise from reprocessing, we might expect the re- 
processing to blur out the soft profile and thus decrease the 
amplitude relative to the hard pulses. The observed ampli- 
tudes might be explained either if the bulk of the emission 
comes from a relatively small reprocessing region (as would 
be suggested by the small inferred 51 = 0.05) and thus limits 
the blurring of the reprocessed profile. Alternatively, if we 
are viewing the hard beam from the neutron star off-axis, but 
the beam strikes the disk directly, then the reprocessing region 
may see pulsations with a significantly higher amplitude than 
what we observe. 

We conclude that the interpretations regarding the disk and 
beam geometries are very different depending on whether we 
consider the £ > 10 keV X-rays which comprise the bulk of 
the surface emission from the neutron star. Further, we note 
that the complex pulse profile and lower pulse fraction in the 
medium band relative to the hard band suggest that the neu- 
tron star emission not be described by a simple pencil beam 
at all energies. Complex beam shapes have been proposed for 
some X-ray pulsars; for example, the hard pulses from Her 
X- 1 have been modeled as a central p encil beam surrounded 
by a fan beam ( iBlum & KrausI 120001) . In our Suzaku obser- 
vations of LMC X-4, the only clear way to explain changes 
in the medium energy pulse profile (particularly the change 



from single to double peaks, and the relatively weak pulsa- 
tions in observation 3) is by intrinsic changes in the beam 
pattern that move the beam out of our line of sight. Weak 
or abse nt pulsations h ave been observed before in the high 
state (La B arbera et alJI200 1). so perhaps the changes in the 
beam pattern are fairly common and happen on timescales of 
^ weeks. 

Despite the difficulty in interpreting the observed pulse pro- 
files, these observations provide further evidence for the soft 
component originating in a reprocessing disk. As mentioned 
in Section [121 using /?bb '- 1.5 x 10** cm and B - 10'^ G, 
we obtain RBB/Rm = 0.7. Keeping in mind that many uncer- 
tainties exist in calculating Rbb and /?„, (such as difficulties 
in measuring B and the poorly understood intera ction of the 
accretion disk with the magnetic field; see, e.g. IScott et al] 
[2000; Hickox et al. 2004), the observed correspondence be- 
tween Rbb and R^ is rather striking and suggests that the hard 
X-rays are reprocessed by optically thick material near the 
magnetosphere. 

Finally, we reiterate that the interpretation of disk re- 
processing for the soft component is pre ferre d on phys- 
ical grounds. As discussed in Section 13.21 a thermal 
bremsstrahlung or MEKAL model (rather than a blackbody) 
adequately describes the soft spectral component in LMC 
X-4, but we can rule out this mechanism because it would 
require too large of an emission region. We conclude, as 
iHickox & Vrtilekl (l2005h have inferred for SMC X-1, the soft 
X-ray component of LMC X-4 is analogous to that of Her X- 
1 and hence due to reprocessing of hard X-rays by the inner 
accretion disk. 

As an aside, we note that if for LMC X-4 we consider the 
"hard" flux only at 2-10 keV, we obtain a smaller Lx and 
thus for Rbb/R,!! ^ 0.3, which is in the range of the values 
{Rbb/Riii ^ 0.2-0.4) measured for Her X -1 over similar en- 
ergy bands (I Scottet al. 2000; Hickox et al. 2004). This sug- 
gests that Rbb may be somewhat closer to /?„, for Her X-1 as 
well, if we consider the full luminosity of the hard component 
which peaks at higher X-ray energies. 

A more detailed understanding of the emission from LMC 
X-4 will require pulse-phase-resolved spectroscopy through- 
out the superorbital period, as well as a more sophisticated 
model of the beam and disk geometry. In the future, the 
International X-ray ObservatorVA with its greater spectral 
resolution and effective area, may enable us to undertake de- 
tailed, high-resolution pulse-phase spectroscopy and extend 
our studies to the low state of the superorbital cycle, enabling 
us to probe further into the accretion mechanism of LMC X-4, 
and gain a better understanding of accretion in highly magne- 
tized environments. 
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